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ABSTRACT 


This  thesis  describes  three  specialized  branch-and-bound  (B&B)  algorithms  for  solving  a 
mixed-integer  program  (MIP)  that  incorporates  standard  “big -M”  constructs.  The  goal  is 
to  identify  valid  values  for  M  that  also  lead  to  short  solution  times. 

One  algorithm  initializes  large  instances  of  M  (giving  a  weak  relaxation  of  the 
MIP),  and  decreases  these  as  required  to  increase  efficiency  of  the  standard  B&B.  Two 
algorithms  initialize  small  and  possibly  invalid  instances  of  M,  and  subsequently  increase 
those  values  in  an  attempt  to  ensure  solution  validity.  Each  algorithm  requires  a  model- 
specific  test  condition  to  detect  weak  or  invalid  M’s. 

We  test  all  algorithms  on  an  uncapacitated  A-median  problem  (a  variant  of  the 
uncapacitated  facility  location  problem),  and  one  algorithm  on  a  shortest-path  interdiction 
problem  (SPIP).  We  observe  substantial  reduction  in  run  times  in  almost  all  cases  tested. 
When  solving  for  exact  solutions,  computational  results  show  that  the  proposed 
algorithms  may  reduce  solution  times  by  up  to  75%  for  the  uncapacitated  A-median 
problem  and  99%  for  the  SPIP.  When  the  algorithms  yield  marginally  suboptimal 
solutions,  substantial  solution-time  improvements  are  also  recorded. 

While  testing  is  limited,  this  thesis  serves  as  a  proof-of-concept  that  the  proposed 
adaptive  algorithms  can  be  effective  in  reducing  solution  times  and  producing  optimal  or 
nearly  optimal  solutions. 
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EXECUTIVE  SUMMARY 


The  continuous  relaxation  of  a  mixed-integer  program  (MIP)  containing  standard  “big- 
M”  constructs  can  be  weak  because  large  values  of  M  (“bound  coefficients”)  are  required 
to  ensure  the  model's  validity.  This,  in  turn,  may  lead  to  slow  convergence  when  solving 
the  model  using  a  standard  branch-and-bound  (B&B)  algorithm.  For  the  purpose  of  this 
thesis,  we  consider  bound  coefficients  to  be  unsuitable  if  they  are  so  large  solution  times 
may  be  long,  or  so  small  that  the  model  may  be  invalid.  One  approach  to  reduce  solution 
times  is  to  pre-calculate  tight  and  valid  bound  coefficients  prior  to  solving  the  MIP.  A 
second  approach,  which  is  the  focus  of  this  thesis,  modifies  the  standard  B&B  algorithm 
to  initialize  bound  coefficients  to  possibly  unsuitable  values,  but  then  attempts  to 
recognize  and  update  those  coefficients  by  either  increasing  or  decreasing  them  as 
required  for  validity  or  short  solution  times,  respectively.  This  process  repeats  until  the 
algorithm  deems  all  bound  coefficients  suitable.  The  algorithms  are  termed  “adaptive” 
because  they  repeatedly  update  bound  coefficients. 

Three  adaptive  algorithms  are  developed.  Two  of  the  algorithms  continually 
update  bound  coefficients  throughout  the  B&B  tree,  whereas  the  other  updates  bound 
coefficients  only  at  the  root  node  of  the  tree,  i.e.,  before  any  branching  occurs.  Briefly, 
the  algorithms  are: 

1.  Algorithm  1.  This  algorithm  initializes  large  (weak)  bound  coefficients 
and  reduces  them  throughout  the  B&B  tree  as  required  for  efficiency. 

2.  Algorithm  2.  This  algorithm  initializes  small  and  potentially  invalid 
bound  coefficients,  but,  at  the  root  node  only,  increases  each  coefficient  that  is 
recognized  by  the  algorithm  to  be  potentially  invalid  before  commencing  the  exploration 
of  the  B&B  tree.  When  updating  occurs,  the  problem  is  re-solved  from  the  root  node 
with  the  updated  bound  coefficients. 

3.  Algorithm  3.  This  algorithm  initializes  small  and  potentially  invalid 
bound  coefficients  but,  throughout  the  B&B  tree,  increases  each  coefficient  that  is 


recognized  by  the  algorithm  to  be  potentially  invalid.  When  updating  occurs,  the 
problem  is  re-solved  from  the  root  node  with  the  updated  coefficients. 

Central  to  all  algorithms  is  the  need  for  a  model-specific  test  condition  that 
signifies  either  potentially  invalid  or  weak  bound  coefficients. 

For  demonstration  purposes,  all  algorithms  are  tested  on  an  instance  of  a  weak 
formulation  of  the  uncapacitated  k-median  problem.  The  problem  comprises  50 
locations,  from  which  5  must  be  chosen  to  fulfill  the  demands  of  100  customers  for  a 
single  product  at  minimum  cost.  Solution  times  for  the  proposed  algorithms  are 
compared  to  solution  times  obtained  when  the  same  formulation  is  solved  using  the 
standard  B&B  with  standard,  provably  valid,  bound  coefficients.  In  the  instances  where 
Algorithms  1  and  2  yield  the  optimal  solution,  we  observe  a  reduction  in  solution  time  of 
up  to  75%  and  62%,  respectively.  Algorithm  3  yields  the  optimal  solution  for  all 
instances  tested,  with  up  to  a  37%  decrease  in  solution  time.  An  alternative  test  condition 
for  Algorithm  3  also  yields  an  optimal  solution  in  all  cases,  but  with  excessive  solution 
times. 

We  test  Algorithm  3  on  an  instance  of  a  shortest-path  interdiction  problem  (SPIP) 
comprising  50  nodes  and  416  arcs.  The  basic  model  seeks  to  maximize  the  shortest-path 
length  in  a  network  by  interdicting  arcs  and  making  them  impassable  to  a  network  user. 
We  solve  a  standard  MIP  formulation  derived  by  taking  the  dual  of  the  inner  shortest- 
path  model.  Algorithm  3  exhibits  up  to  a  99.8%  reduction  in  solution  time.  We  also 
observe  a  maximum  3%  deviation  from  the  optimal  objective  value. 

This  thesis  serves  as  a  proof-of-concept  to  show  that  the  proposed  adaptive 
algorithms  can  be  used  to  substantially  reduce  solution  times  of  MIPs  containing  standard 
big-M  constructs.  Note  that,  since  the  solution  quality  of  all  proposed  algorithms  cannot 
be  guaranteed  a  priori  for  all  problem  instances,  these  adaptive  algorithms  are  heuristics. 
Although  results  are  promising,  further  research  must  be  conducted  before  practical 
implementation  can  be  adopted. 
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I.  INTRODUCTION 


This  thesis  describes  specialized  branch-and-bound  (B&B)  algorithms  for  solving 
a  mixed-integer  program  (MIP)  in  which  binary  variables  x;  control  continuous  variables 

yt  representing  “activity  levels”:  For  each  i,  either  both  variables  have  values  0,  or  the 
controlling  binary  variable  becomes  1  and  the  continuous  activity  is  then  allowed  to  take 
on  positive  values,  up  to  the  value  of  its  bound  coefficient.  In  particular,  we  are 
interested  in  MIPs  that  contain  constraints  of  the  form: 


y-M 

< 

o, 

Vie/ 

(1.1) 

T; 

> 

o, 

Vie/ 

(1.2) 

V 

e 

{0 

,1},  Vie/, 

(1.3) 

where  Mt>  0 ,  for  all  i  el .  The  model  coefficient  Mi  is  the  bound  coefficient  for  the 
continuous  activity  level  yt ,  and  x  is  the  controlling  binary  variable.  Note  that  y,  can 
actually  correspond  to  an  affine  function  of  other  continuous  variables.  The  key 
characteristic  of  the  choice  of  Mi  is  that  this  bound  coefficient  must  be  chosen  such  that 

(1.1)  is  redundant  when  xt  is  1. 

Examples  of  models  that  may  contain  “big -M"  constructs  are  weak  formulations 
of  the  uncapacitated  facility  location  problem  (UFLP),  as  described  by  Khumawala 
(1972),  and  the  closely  related  uncapacitated  k-median  problem  (Lin  &  Vitter,  1992). 
Similar  constructs  appear  in  production-scheduling  problems,  e.g.,  Eppen  and  Martin 
(1987),  Maes  and  Van  Wassenhove  (1988),  and  MIP  fonnulations  of  min-max  (or  max- 
min)  models  such  as  network-interdiction  problems,  e.g.,  Brown  et  al.  (2006),  Wood 
(1993).  In  all  problems,  “variable  upper  bounding  constraints,”  i.e.,  equation  (1.1), 
connect  binary  and  continuous  variables. 

The  continuous  relaxation  of  a  “big -M  model”  can  be  weak  because  large  values 
of  M  are  required  to  ensure  the  model’s  validity.  Some  textbooks  such  as  Markland  and 
Sweigart  (1987,  p.  505)  advocate  using  “a  large  number”  for  bound  coefficients  in  the 
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model.  However,  Camm  et  al.  (1990)  note  that,  by  making  bound  coefficients  large,  the 
feasible  region  of  the  continuous  relaxation  is  expanded  unnecessarily.  This,  in  turn,  may 
lead  to  slow  convergence  when  solving  the  model  using  a  standard  B&B  algorithm. 
Other  textbooks,  such  as  Rardin  (1998,  p.  642),  require  bound  coefficients  to  be 
“sufficiently  large,”  i.e.,  not  excessively  large  for  efficiency,  but  large  enough  for 
validity.  Rardin  emphasizes  that  “whenever  a  discrete  model  requires  sufficiently  large 
big  M’s,  the  strongest  relaxations  will  result  from  models  employing  the  smallest  valid 
choice  of  those  constants.”  For  the  purpose  of  this  thesis,  bound  coefficients  are 
considered  unsuitable  if  they  are  so  small  that  the  model  may  not  be  valid,  or  so  large  that 
solution  times  are  excessive. 

One  approach  to  reduce  solution  times  for  a  “big-M  model”  is  to  predefine  valid 
values  for  the  bound  coefficients  prior  to  solution,  while  attempting  to  make  those  values 
as  tight  as  possible.  Crowder  et  al.  (1983),  Van  Roy  and  Wolsey  (1987)  and  Johnson  et 
al.  (1985)  provide  preprocessing  techniques  to  develop  tighter  values  for  bound 
coefficients.  These  early  techniques  primarily  involve  careful  model  formulation,  and  are 
widely  considered  as  standard  present-day  modeling  techniques.  For  example,  in  an 
UFLP,  using  the  sum  of  demands  for  bound  coefficients  is  more  efficient  than  using  an 
arbitrarily  large  value.  These  tighter  fonnulations,  however,  may  still  be  inefficient  for 
solving  large  problems.  We  note  that  there  may  be  stronger  preprocessing  techniques 
available  that  are  unknown  to  this  author. 

This  thesis  explores  a  second  approach  to  reduce  solution  times  for  a  “big-M 
model.”  The  goal  is  to  modify  the  standard  B&B  algorithm  to  initialize  bound 
coefficients  to  possibly  unsuitable  values,  but  then  to  recognize  and  update  those 
coefficients  by  either  increasing  or  decreasing  them  as  required  for  validity  or  efficient 
solution,  respectively.  We  propose  three  adaptive  algorithms;  two  algorithms  continually 
update  bound  coefficients  throughout  the  B&B  tree,  whereas  the  other  updates  bound 
coefficients  only  at  the  root  node  of  the  tree,  i.e.,  before  any  branching  occurs.  The 
algorithms  are  termed  “adaptive”  because  they  repeatedly  update  bound  coefficients.  The 
solution  quality  of  the  algorithms,  at  the  time  of  writing,  cannot  be  guaranteed  a  priori. 
Hence,  these  adaptive  algorithms  are  heuristics. 
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A.  STANDARD  BRANCH-AND-BOUND 

B&B  implicitly  searches  the  entire  feasible  region  for  the  continuous  relaxation  of 
the  MIP  it  is  solving.  Every  branch  partitions  the  locally  restricted  feasible  region  into 
three  parts,  one  of  which  cannot  contain  an  integer  solution  and  is  therefore  not  explored. 
Assuming  a  minimization  problem,  every  existing  partitioned  region  with  an  objective 
value  that  exceeds  the  objective  value  of  any  known  integer  solution  cannot  contain  the 
optimal  solution;  therefore  it  is  terminated  or  fathomed.  Similarly,  partitioned  regions 
that  are  infeasible  are  also  fathomed.  A  lower  bound  is  also  derived  from  the  smallest 
objective  value  of  continuous  solutions  from  the  existing  partitioned  regions.  The 
partitioning  continues  until  an  integer  solution  is  found  such  that  its  objective  value  is 
equal  to  the  lower  bound.  See  Lawler  and  Wood  (1966)  for  a  survey  of  the  B&B 
algorithm  and  its  applications. 

B.  LITERATURE  REVIEW 

The  adaptive  algorithms  proposed  in  this  thesis,  while  not  exactly  analogous,  are 
similar  in  flavor  to  variable -bound  tightening  techniques  used  in  multi-start  methods  for 
solving  non-linear  programs  (NLPs);  see  Chinneck  (2002).  Multi-start  methods  for  an 
NLP  aim  to  provide  an  optimization  algorithm  with  an  initial  solution  from  which  it  can 
reach  feasibility  efficiently.  These  methods  furnish  random  starting  solutions  within  a 
region  provided  by  the  variable-bounds.  However,  the  variable-bounds  provided  by  the 
original  model  may  define  a  region  that  is  much  larger  than  required  given  the  other 
constraints  in  the  problem.  Therefore,  it  is  desirable  to  define  a  region  that  is  not 
excessively  larger  than  the  region  of  interest.  Tightening  of  variable-bounds  is  then 
applied  to  give  the  multi-start  algorithms  a  smaller  region  from  which  a  good  initial 
solution  is  likely  to  be  found  efficiently.  Chinneck  notes  that  the  probability  of  obtaining 
the  tightest  possible  variable-bounds  that  do  not  eliminate  feasible  solutions  is  small.  He 
proposes  bound-tightening  techniques  via  the  following  means: 
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1.  Depending  on  the  modeler’s  aims  and  analysis  of  the  problem,  the  bounds 
may  be  adjusted  manually  to  contain  suspected  optimal  points  or  simply  to  exclude 
regions  that  are  unlikely  to  contain  an  optimal  solution. 

2.  Presolving  techniques,  such  as  those  proposed  by  Brearley  et  al.  (1975), 
can  be  used  to  tighten  the  variable-bounds  or  even  fix  the  value  of  certain  variables. 

3.  Linear  and  non-linear  interval  analysis  can  be  applied  when  the  bounds 
can  be  inferred  from  constraints  present  in  the  model.  The  idea  is  to  determine  the 
maximum  possible  range  for  each  variable  v,  by  considering  the  constraints  where  v 
appears  and  the  bounds  on  other  variables  present  in  those  constraints.  This  is  also  a 
useful  method  to  detect  infeasibility  of  a  model  should  the  maximum  possible  range  for 
any  variable  fall  outside  the  stated  range  in  the  original  problem. 

Chinneck  also  proposes  an  alternative  method  to  bound  variables  in  an  NLP.  This 
method  starts  with  a  small  region  and  expands  it  in  several  stages.  However,  it  is 
possible  that  the  final  bounds  derived  may  be  invalid,  resulting  in  suboptimal  solutions. 

Chinneck  describes  how  variable-bounds  can  be  increased  and  decreased  to  obtain 
a  good  initial  point  efficiently.  In  a  similar  fashion,  an  adaptive  B&B  algorithm  will 
increase  or  decrease  bound  coefficients  to  help  improve  efficiency;  rules  for  changing 
coefficients  will  attempt  to  ensure  validity. 

C.  ADAPTIVE  BRANCH  AND  BOUND 

The  adaptive  B&B  algorithms  described  in  this  thesis  explore  these  three  distinct 
methods: 

1 .  Initialize  each  bound  coefficient  to  a  large  and  probably  weak  value  pre¬ 
calculated  via  standard  techniques.  Such  coefficients  are  henceforth  termed  baseline 
coefficients.  Then,  during  the  exploration  of  the  B&B  tree,  reduce  each  coefficient  that  is 
recognized  by  the  algorithm  to  be  potentially  weak. 

2.  Initialize  each  bound  coefficient  to  a  small  value  that  is  potentially  invalid. 
Then,  at  the  root  node  only,  increase  each  coefficient  that  is  recognized  by  the  algorithm 

to  be  potentially  invalid  before  commencing  the  exploration  of  the  B&B  tree. 
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3.  Initialize  each  bound  coefficient  to  a  small  value  that  is  potentially  invalid. 
Then,  throughout  the  B&B  search,  increase  each  coefficient  that  is  recognized  by  the 
algorithm  to  be  potentially  invalid.  Whenever  updating  occurs,  the  problem  is  re-solved 
from  the  root  node  with  the  updated  bound  coefficients. 

The  algorithms  recognize  a  potentially  unsuitable  bound  coefficient  at  a  node 
when  a  test  condition  is  met,  and  attempt  to  increase  or  decrease  the  coefficient 
depending  on  the  method  being  implemented.  Test  conditions  are  model-specific,  and 
when  met,  signify  that  a  coefficient  is  potentially  invalid  or  weak.  A  test  condition  that 
identifies  a  potentially  invalid  bound  coefficient  is  referred  to  as  an  “invalidity  test 
condition,”  and  a  test  condition  that  identifies  a  potentially  weak  bound  coefficient  is 
referred  to  as  an  “inefficiency  test  condition.”  Once  the  potentially  unsuitable  bound 
coefficients  are  updated,  the  adaptive  algorithm  continues  until  a  set  of  bound 
coefficients  is  deemed  unsuitable  at  a  subsequent  step  in  the  algorithm. 

The  adaptive  B&B  algorithms  are  described  in  detail  in  Chapter  II,  and  are 
applied  to  a  weak  formulation  of  the  uncapacitated  k-median  problem  in  Chapter  III  and  a 
shortest-path  interdiction  problem  (SPIP)  in  Chapter  IV.  (The  aforementioned  weak 
formulation  would  probably  not  be  used  in  practice,  and  is  used  in  this  thesis  only 
because  it  provides  a  simple,  paradigmatic  test  problem.)  This  thesis  conducts  a  proof- 
of-concept  study  to  show  that  the  proposed  adaptive  algorithms  can  be  used  to  reduce 
solution  times  of  MIPs  containing  constraints  of  the  form  (1.1)  -  (1.3).  Positive  results 
on  the  test  problems  will  indicate  that  further  research  in  this  area  is  warranted. 


5 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


6 


II.  ADAPTIVE  BRANCH-AND-BOUND  ALGORITHMS 


This  chapter  describes  the  three  variants  of  adaptive  B&B  algorithms  in  further 

detail. 

The  following  notation  is  common  to  all  three  algorithms: 
n  node  in  the  B&B  tree;  n  -  0  is  the  root  node 
M\  initial  value  of  bound  coefficient  used  in  the  /-th  constraint 

Mi  adjusted  value  of  bound  coefficient  used  in  the  i- th  constraint 

Vn  index  set  of  “violated”  constraints,  i.e.,  constraints  whose  coefficients  Mi 
meet  the  inefficiency  or  invalidity  test  condition,  at  node  n. 

Note  that  the  details  of  the  inefficiency  and  invalidity  test  conditions  are  not 
specified  until  the  test  models  are  introduced  in  Chapters  III  and  IV. 

A.  ALGORITHM  1:  REDUCING  A  LARGE,  INITIAL  M 

Algorithm  1  assumes  that  large,  valid  baseline  bound  coefficients  have  been  pre¬ 
calculated  via  standard  techniques.  The  algorithm  seeks  to  recognize  weak  coefficients 
and  decrease  their  values  appropriately.  The  process  of  updating  bound  coefficients  can 
occur  anywhere  in  the  B&B  tree. 

We  define  parameter  /?,  0  </?<!,  as  the  coefficient  reduction  factor,  where 
1  -  f  is  the  factor  by  which  a  bound  coefficient  is  reduced  each  time  an  inefficiency  test 
condition  is  met.  Algorithm  1  follows: 
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Algorithm  1 

Input: 

A  MIP  containing  constraints  of  the  form  (1.1),  (1.2)  and  (1.3). 

Coefficient  reduction  factor  (5  ,  0  <  (3  <  1 . 

Inefficiency  test  condition  with  associated  parameters.  (See  equation  3.B.1) 

Baseline  coefficients  M\,  V/  e  /  . 

Output: 

A  feasible  solution  to  the  MIP. 

Step  0 

Start  a  standard  B&B  algorithm; 

Solve  root  node  n  =  0 ; 

Step  1 

Select  an  active  node  n  to  process  using  the  B&B  algorithm’s  selection 

procedure; 

If  node  n  is  fathomed  then  go  to  step  4; 

Step  2 

Determine  I'n  via  inefficiency  test  condition; 

If  Pn  =  0  then  let  the  B&B  algorithm  continue,  possibly  creating  child 

nodes,  and  then  proceed  to  step  1 ; 

If  I'n  *  0  then  go  to  step  3; 

Step  3 

Let  M.  <—  f3deptH,,)Ml,  V/  g  F ,  where  depth(n)  is  the  depth  of  the  node  n 

being  solved  within  the  B&B  tree; 

Impose  Mn  V/  el  on  any  descendants  of  node  n; 

Step  4 

If  there  exists  an  active  node  in  the  B&B  tree  then  go  to  step  1 ; 

Step  5 

Print  solution; 

Stop; 

End  of  Algorithm  1 

Remark:  By  observation,  a  larger  value  for  /?  implies  that  M j  is  reduced  at  a 


slower  rate  than  if  a  smaller  / 3  is  used. 
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Ideally,  the  incumbent  value  for  the  bound  coefficient  would  be  used  in  the 
updating  process,  i.e.,  Mj  <—  j3Mi .  However,  this  value  could  not  be  retrieved  during  the 

execution  of  the  B&B  algorithm  using  Xpress  software  (see  Section  D).  Hence,  scaling 
according  to  depth  in  the  B&B  tree  serves  as  a  proxy. 

The  reader  should  note  that  the  algorithm  is  a  heuristic  because  any  bound 
coefficient  may  be  reduced  below  its  smallest  valid  value. 

B.  ALGORITHM  2:  INCREASING  A  SMALL,  INITIAL  M  AT  ROOT  NODE 

The  second  algorithm  attempts  to  identify  invalid  bound  coefficients  based  on  the 
root-node  solution,  that  is,  before  any  branching  occurs.  Using  a  small  and  possibly 
invalid  initial  value  for  each  bound  coefficient,  the  relaxed  problem  is  solved,  and  if  the 
algorithm  detects  potentially  invalid  coefficients,  these  coefficients  are  increased  and  the 
root  node  is  re-solved.  This  process  continues  until  the  algorithm  does  not  detect  any 
potentially  invalid  bound  coefficients.  When  this  occurs,  the  standard  B&B  algorithm 
proceeds  normally  (without  further  updates  of  the  coefficients). 

We  define  y>\  to  be  the  increase  factor  of  the  bound  coefficient,  and  update 
Mt  <-  yMj  for  all  i  e  70v .  Algorithm  2  follows: 
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Algorithm  2 

Input: 

A  MIP  containing  constraints  of  the  form  (1.1),  (1.2)  and  (1.3). 

Increase  factor  y>  1  (for  processing  the  root  node  n  —  0  ). 

Invalidity  test  condition  with  associated  parameters  (for  processing  the  root 

node/7  =  0).  (See  equation  3. B. 2) 

Small  value  of  M\ ,  Vie/. 

Output: 

A  feasible  solution  to  the  MIP  defined  with  weak  bound  coefficients. 

Step  0 

M,  <-  Mf,  Vizi; 

Step  1 

Solve  root  node  n  =  0  ; 

Step  2 

Determine  /„  via  invalidity  test  condition; 

If  I'0  =  0  then  go  to  step  4; 

If  Iv0  *  0  then  go  to  step  3; 

Step  3 

Set  Mt  <-  yMi ,  V?e/0V; 

Go  to  step  1 ; 

Step  4 

Implement  a  standard  B&B  algorithm; 

Step  5 

Print  solution; 

Stop; 

End  of  Algorithm  2 

This  algorithm  is  optimistic  as  it  relies  only  on  root  node  solutions  in  the  attempt 


to  derive  valid  bound  coefficients. 

C.  ALGORITHM  3:  INCREASING  A  SMALL,  INITIAL  M  WITHIN  THE 
SEARCH  TREE 

In  Algorithm  2,  bound  coefficients  are  increased  at  the  root  node  only.  The  last 
algorithm.  Algorithm  3,  allows  bound  coefficients  to  be  increased  anywhere  in  the  B&B 
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search  tree,  including  the  root  node.  To  do  this.  Algorithm  3  leverages  on  Algorithm  2  to 
process  the  root  node  n  =  0 ;  other  inputs  are  defined  as  required  to  process  tree  nodes 
n  >  1 .  Should  Algorithm  3  detect  that  incumbent  bound  coefficients  are  possibly  invalid, 
these  coefficients  are  increased  and  the  problem  is  re-solved  from  the  root  node.  This 
process  continues  until  the  standard  B&B  algorithm  converges  with  all  bound 
coefficients  deemed  valid. 

Define  (p>  1  to  be  the  increase  factor  in  the  B&B  tree  ( n  >  1 ),  similar  to  how  y 
is  defined  in  Algorithm  2.  Algorithm  3  follows: 
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Algorithm  3 

Input: 

Same  inputs  as  for  Algorithm  2  (for  processing  the  root  node  n-  0  ). 

Increase  factor  (p>\  (for  processing  tree  nodes  n  >  1). 

Invalidity  test  condition  with  associated  parameters  (for  processing  tree  nodes 

n>\)  (See  equations  3.B.3,  3.B.4  and  4.B.1) 

Output: 

A  feasible  solution  to  the  MIP  defined  with  weak  bound  coefficients. 

Step  0 

Implement  Algorithm  2  until  70v  =  0  ,  i.e.,  until  standard  B&B  proceeds 

normally; 

Step  1 

Select  an  active  node  n  to  process  using  the  B&B  algorithm’s  selection 

procedure; 

Determine  Pn  via  invalidity  test  condition; 

If  f'n  =  0  then  go  to  step  3; 

If  Pn  *  0  then  go  to  step  2; 

Step  2 

Set  M i  <—  <pM , ,  Vie/;; 

Go  to  step  0; 

Step  3 

If  there  exists  an  active  node  in  the  B&B  tree  then  go  to  step  1 ; 

Step  4 

Print  solution; 

Stop; 

End  of  Algorithm  3 

As  opposed  to  Algorithm  2,  which  operates  only  at  the  root  node,  this  algorithm 
provides  the  flexibility  to  adjust  bound  coefficients  as  the  search  tree  grows.  This 
increases  the  likelihood  of  attaining  an  optimal  solution  as  the  bound  coefficient  values 
required  for  validity  might  not  be  obtained  by  iterating  at  the  root  node  alone. 
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D.  IMPLEMENTATION  OF  ALGORITHMS 


All  models  and  algorithms  are  implemented  using  Xpress-IVE  Version  1.22.04  as 
the  solving  engine.  Computations  are  performed  on  an  Acer  Aspire  474 1G  laptop  with  6 
gigabytes  of  RAM,  running  on  an  Intel1'  Core™  i5-M460  processor  (2.53GHz,  3MB  L3 
cache)  and  Windows8  7  operating  system. 

For  the  purpose  of  evaluating  the  performance  of  the  adaptive  algorithms  against 
a  standard  B&B  algorithm,  we  disable  all  of  Xpress-IVE ’s  presolve  techniques,  cutting- 
planes  strategies,  and  heuristics. 
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III.  TEST  MODEL  1:  UNCAPACITATED  ^-MEDIAN  PROBLEM 


This  chapter  describes  the  weak  formulation  of  the  uncapacitated  A-median 
problem  and  solves  it  using  Algorithms  1,  2  and  3.  We  state  the  test  conditions  for  and 
analyze  the  results  of  each  algorithm. 

A.  MODEL  FORMULATION 

In  the  uncapacitated  A-median  problem,  facilities  of  unrestricted  size  are  placed 
among  n  possible  locations  with  the  objective  of  minimizing  cost  for  satisfying  given 
demands  for  a  single  product  at  m  customer  locations.  Costs  include  a  per-unit  shipping 
cost  for  product  shipped  from  facility  i  to  customer  j,  and  a  per-unit  penalty  at  each 
customer  location  for  any  unmet  demand.  The  number  of  facilities  to  be  opened  is  given. 

We  note  that  the  uncapacitated  A-median  problem  is  a  variant  of  the  UFLP.  The 
difference  is  that  UFLP  applies  a  fixed  cost  for  opening  any  facility,  rather  than  limiting 
the  number  of  opened  facilities.  See  Reese  (2005)  for  a  summary  of  solution  methods  for 
the  uncapacitated  A-median  problem. 

As  highlighted  in  Chapter  I,  the  uncapacitated  A-median  problem  formulation 
presented  here  is  a  weak  version  that  is  used  for  testing  purposes;  this  formulation  might 
not  be  used  in  practice.  (Efroymson  &  Ray  1966  do  study  the  weak  formulation, 
however.)  The  differences  between  the  weak  and  strong  fonnulations  are  discussed  at  the 
end  of  this  section. 

The  mathematical-programming  fonnulation  of  the  uncapacitated  A-median 
problem  follows: 

Indices  and  Index  Sets 
i  <e  /  facilities  (/  =  {*),..., in }  ) 

J  e  J  customers  (  J  =  { j\ ,  ■  ■  ■ ,  jm  }  ) 
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Parameters 


d j  units  of  demand  at  customer  j  for  some  commodity 

M i  bound  coefficient  for  total  shipping  from  facility  i 

Cy  per-unit  shipping  cost  from  facility  i  to  customer  j 

Pj  per-unit  penalty  of  unmet  demand  at  customer  j 
k  exact  number  of  facilities  that  must  be  opened 

Remark:  At  convergence  of  the  adaptive  algorithms,  the  bound  coefficients  M  i  should 

be  large  enough  such  that  no  unmet  demand  arises.  However,  we  acknowledge  that  some 
of  the  bound  coefficients  may  be  invalid,  at  least  temporarily,  as  the  adaptive  B&B 
algorithms  proceed. 

Variables 

x  1  if  facility  i  is  opened,  0  otherwise 

ytj  number  of  units  shipped  from  facility  i  to  customer  j 


Sj  units  of  unmet  demand  at  customer  j 
Formulation 

2*  =  1™  Ilvrl/'r''/  (3.A.1) 

i  j  j 

s,t'  X-fy  <  Mtxt,  V/ el  (3.A.2) 

j 

Z Fy  +  SJ  >  dr  VjeJ  (3.A.3) 

i  J 

Xv  =  k  (3.A.4) 

i 

y  >  0  (3.A.5) 
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s  >  0 


(3.A.6) 


x  e  {0,1}"*”’  (3.A.7) 

The  objective  function  (3.A.1)  defines  the  total  cost  to  be  minimized.  This 
comprises  transportation  costs  and  penalties  should  any  unfulfilled  demand  occur. 
Constraints  (3. A. 2)  govern  the  relationship  between  production  and  setup  of  a  facility. 

Constraints  (3 .A. 3)  determine  the  amount  of  unfulfilled  demand  and  charge  the 
corresponding  penalty  to  the  objective  function  (3.A.1).  Note  that,  while  fixing  the 
number  of  facilities  to  be  opened  in  (3. A. 4)  makes  this  problem  an  uncapacitated  k- 
median  problem,  we  believe  the  proposed  algorithms  can  be  applied  to  weak  UFLP 
formulations,  also. 

Remark:  The  strong  formulation  presented  by  Balinski  (1966)  replaces 

constraints  in  (3. A. 2)  with  y Id /  <  x.  for  all  i  el,  and  j  &J ,  where  y  Id .  represents  the 
fraction  of  customer  f  s  demand  that  is  satisfied  by  facility  /.. 

B.  TEST  CONDITIONS  FOR  UNCAPACITATED  ^-MEDIAN  PROBLEM 

This  section  discusses  possible  test  conditions,  specific  for  the  uncapacitated  k- 
median  problem,  to  detect  potentially  weak  or  invalid  bounds.  We  first  present  data  for  a 
small  problem  that  is  used  to  explain  the  implementation  of  the  algorithms.  Actual 
computational  tests  are  carried  out  on  a  larger  problem  (see  Section  C). 

1.  Illustrative  Example  Data 

Consider  four  possible  facility  locations  and  five  customers.  Two  facilities  must 
be  opened.  Each  of  the  five  customers’  demands  is  shown  in  Table  1. 


Customer 

1 

2 

3 

4 

5 

Demand,  d . 

180 

240 

210 

300 

150 

Table  1.  Customer  demands 
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Table  2  presents  per-unit  transportation  costs  ctj  for  shipments  from  facility  i  to 
customer  j. 


'x\sCustomer 

2 

3 

4 

5 

F  acility^^^^ 

1 

2 

5 

3 

2 

2 

2 

1 

3 

5 

4 

2 

3 

5 

3 

1 

4 

3 

4 

2 

5 

2 

1 

3 

Table  2.  Per-unit  transportation  costs 


For  each  customer,  penalty  costs  per  unit  of  unmet  demand  are  set  to  be  higher 
than  the  maximum  of  all  transportation  costs.  Specifically,  we  set  p .  -  6  for  all  j. 

For  this  problem,  an  optimal  solution  is  x*  =  (0,  1,  0,  1)  with  an  optimal  objective 
value  z*  =  1920. 

2.  Algorithms  with  Test  Conditions 

a.  Test  Condition  for  Algorithm  1 

Starting  with  the  baseline  coefficient,  Mt  <—  M1 ,  we  infer  that  a  bound 
coefficient  is  possibly  weak  whenever  the  total  demand  satisfied  by  facility  i  is  less  than  a 
predefined  proportion  of  M. .  Let  depth(n)  denote  the  depth  of  node  n  in  the  B&B  tree 

being  processed.  The  inefficiency  test  condition  for  Algorithm  1  applied  to  the 
uncapacitated  k-median  problem  is 

i  e  /;  if  Z  J  y0  <  aj3dep,Hn)~lMl ,  (3  -B • 1 ) 

where  a  <  1  is  the  inefficiency  test  factor  that  defines  the  aforementioned  proportion. 
After  applying  the  test  condition  to  compute  /j ,  we  let  M.  <—  Jflepth{") ,  for  all  /  e  Vn , 
and  Algorithm  1  continues. 
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We  note  that,  if  the  incumbent  value  for  the  bound  coefficient  could  be 
retrieved  in  Xpress-IVE,  the  test  condition  above  would  be  replaced  by  ^  .  yt]  <  aMi  , 
where  M;.  is  the  relevant  bound  coefficient.  We  also  note  that,  in  order  for  the  algorithm 
to  reduce  a  bound  coefficient  Mi ,  the  total  flow  from  facility  i  cannot  exceed  M,  a . 
Hence,  a  small  inefficiency  test  factor  a  implies  a  smaller  Vn  and  is  a  restriction 
compared  to  a  larger  value  of  a  .  The  larger  ft  is,  the  slower  the  coefficient  is  reduced. 

To  illustrate,  we  now  solve  the  uncapacitated  ^-median  example  using 
Algorithm  1  with  the  following  inputs: 

Input:  Uncapacitated  ^-median  problem  defined  in  previous  section. 

Coefficient  reduction  factor  /3  =  0.9 . 

M'=Tjjdj  =  1080,  Vie/ 

Inefficiency  test  condition  X  ,  ^  aj3depth(n)~]  M ' ,  with  a  =  0.2 

Step  0:  A  standard  B&B  algorithm  solves  the  root  node.  The 

solution  is  xr=(l.00,  0.528,0.194,0.278),  and  the  amount  each  facility  supplies  to 
customers  is  shown  in  Table  3. 

Step  1 :  Since  only  the  root  node  has  been  solved  and  the  solution  is 

continuous,  it  is  not  fathomed.  Algorithm  1  proceeds  to  process  the  root  node. 

Step  2:  Facilities  i  =  1  and  i  —  3  meet  the  inefficiency  test 

condition  (see  Table  3),  thus  Vo  =  {1,3} . 

Step  3:  Let  Mi  <-  fiMj  =  0.9  x  1080  -  972  ,  Vz  e  70v ,  and  impose 

the  updated  values  on  all  descendant  nodes  in  the  B&B  tree. 

Step  4:  The  B&B  algorithm’s  selection  procedure  selects  one  of  the 

two  active  nodes  (i.e.,  children  nodes  of  the  root  node)  to  process,  and  the  adaptive 
algorithm  continues. 
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z*= 1920. 


Continuing  this  process,  an  optimal  solution  obtained  is  x*  =  (0,  1,  0,  1), 


^\Qustomer 

Facility  ^  v 

1 

2 

3 

4 

5 

Total 

I.n.<216? 

1 

150 

Yes 

2 

180 

240 

No 

3 

210 

Yes 

4 

300 

No 

Table  3.  Root  node  solution  for  flow  variables  (Algorithm  1) 


b.  Test  Condition  for  Algorithm  2 

When  starting  with  a  small  bound  coefficient,  Algorithm  2  detects 
potential  invalidity  whenever  the  solution  at  the  root  node  shows  that  the  total  demand 
satisfied  by  facility  i  is  at  least  a  predefined  proportion  8  of  the  incumbent  bound 
coefficient.  The  invalidity  test  condition  for  Algorithm  2  applied  to  the  uncapacitated  k- 
median  problem  is 

(3.B.2) 

where  8  <  1  is  the  invalidity  test  factor.  After  applying  the  test  condition  to  compute  I'0 , 
we  let  M.  <—  yM t  ,  for  all  i  e  70v ,  and  the  problem  is  restarted.  Furthermore,  we  set 

Mt  <—  X  j  dj ,  for  all  i  e  70v  such  that  yM .  >  X ,  d . .  (In  a  similar  argument  used  for  the 
test  factor  a  in  Algorithm  1,  a  larger  8  here  implies  a  more  stringent  condition  to  signify 
possible  invalidity  of  M.  .) 

To  illustrate,  we  now  solve  the  uncapacitated  ^-median  example  using 
Algorithm  2  with  the  following  inputs: 
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Input:  Uncapacitated  ^-median  problem  defined  in  previous  section. 

Increase  factor  y  =  0.1  applied  at  root  node. 

Invalidity  test  condition  ytj  ^  SM  t  ,  with  5  -  0.3 . 

M'  =Tjjdj/k  =  540,  Vie/. 

The  rationale  for  using  the  average  demand  fulfilled  by  an  opened  facility  as  the  initial 
value  for  bound  coefficients  is  discussed  in  a  subsequent  section. 

Step  0:  Set  M.  M\,  for  all  i  e  7. 

Step  1:  The  solution  to  the  root  node  is 

xr  =(0.278,  0.778,  0.389,  0.556)  and  the  amount  each  facility  supplies  to  the  customers 
is  shown  in  Table  4. 

Step  2:  Facilities  i  =  2,  3,  4  meet  the  invalidity  test  condition  (see 

Table  4),  and  thus  70v  =  {2, 3, 4}  . 

Step  3:  Let  Mt  <-  yMj  =1.1x540  =  594,  Vi  e  IVQ  . 

The  root  node  is  solved  again  using  the  updated  values  for  M2 ,  M3  and 
M4 ,  and  checked  using  the  invalidity  test  conditions.  Continuing  this  process,  an  optimal 
solution  obtained  is  x*  =  (0,  1,  0,  1),  z*=1920. 


\Customer 

2 

3 

4 

5 

Total 

L,4„>162? 

Facility  \ 

1 

150 

150 

No 

2 

180 

240 

420 

Yes 

3 

210 

210 

Yes 

4 

300 

300 

Yes 

Table  4.  Root  node  flow  variables  solution  to  MIP  defined  with  invalid  bound 
coefficients  (Algorithm  2) 
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c.  Test  Condition  for  Algorithm  3 

The  invalidity  test  condition  for  Algorithm  3,  applied  to  the  uncapacitated 
/i-media  problem,  is  similar  to  that  of  Algorithm  2: 

i e  K  if  Hj-y.j  ^  (!)M,  >  (3-B-3) 

where  (j)<  1  is  the  invalidity  test  factor  (similar  to  5  in  Algorithm  2’s  invalidity  test 
condition).  As  we  have  described  how  Algorithm  2  processes  node  n-  0 ,  we  proceed  to 
describe  the  processing  of  nodes  n>  1 .  After  applying  the  test  condition  to  compute  I'n , 

we  let  M,  <—  (pMi  for  all  i  e  Fn  and  the  problem  is  restarted  at  n  =  0 .  Furthermore,  we 

set  M <—  Zj  d . ,  for  all  i  e  F  ,  such  that  M  ■>2_,  .d  ■ . 

Using  y  =  ^  =  1.2,  <5  =  0. 1  and  ^  =  0.9,  the  example  is  solved  using 

Algorithm  3.  An  optimal  solution  obtained  is  x*  =  (0,  1,  0,  1),  z*=T920. 

d.  Alternative  Test  Condition  for  Algorithm  3 

We  now  present  an  alternative  method  for  detecting  potential  invalidity  in 
Algorithm  3.  We  showed  earlier  how  bound  coefficients  in  this  algorithm  are  deemed 
potentially  invalid  when  the  total  demand  satisfied  by  facility  i  is  at  least  a  predefined 
proportion  of  the  bound  coefficient.  A  second  method  of  detecting  possible  infeasibility 
of  bound  coefficients  lies  in  the  unique  structure  of  an  uncapacitated  ^-median  problem: 
Assuming  no  two  transportation  costs  to  a  given  customer  are  identical,  i.e.,  ctj  =£  c(j ,  for 

all  i,  i',  i  e  /,  i  ^  i' ,  for  all  j  e  J  ,  then  at  any  node,  if  the  total  number  of  positive  flows 
from  facilities  to  customers  is  not  equal  to  the  number  of  customers,  then  the  bound 
coefficients  are  invalid.  This  condition  implies  that  that,  for  any  instance  of  the  problem 
without  two  identical  costs,  an  optimal  solution  to  every  node  requires  each  customer  to 
be  serviced  by  exactly  one  facility. 
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Lemma  1:  At  each  node,  let  S  be  the  subset  of  J  such  that  j  e  S  if  and  only  if  y*  >  0, 
where  y*  is  the  vector  of  optimal  flow  from  facilities  to  customers.  If  ctJ  *  cfj ,  for  any 
i*i'  el ,  and  any  j  e  J  ,  then  |5”|  =  | j\ . 


Proof.  By  contradiction,  suppose  |^|  >  |,/| .  Then  there  exists  at  least  one  customer  whose 
demand  is  satisfied  by  at  least  two  facilities,  i.e.,  3  j  such  that  {y*~  \  y*~  >  0}  >  2 . 
Because  ctj  ^  c fJ  for  any  i  *  i'  s  /,  and  any  j  &  J  ,  we  can  reduce  the  objective  value  by 
increasing  the  bound  coefficient  of  facility  i  where  /  =  argmin  \c  |  y*~.  >  0J ,  and 


assigning  all  demand  d~  to  facility  i .  Therefore,  y"  is  not  optimal. 


Now,  suppose  \S\  <  |  J\ .  This  implies  that  the  total  capacity  of  the  opened  facilities  is  less 

than  the  total  demand,  and  hence  unmet  demand  occurs  for  some  customers.  Since 
p.  >cij  for  all  iel,  j  e  J ,  the  objective  value  can  be  reduced  by  increasing  bound 

coefficients.  Therefore,  y*  is  not  optimal. 

Hence,  at  optimality,  \S\  -  \j\  holds.  Q.E.D. 

Remark  1:  If  there  exists  i  ^  i'  with  c.  -c,~  for  some  /  g  J ,  and 

ij  i  j  J 

{/,/'}  cz  argminjc.,-.  |  y*~.  >  0},  then  we  can  let  the  demand  at  customer  j  be  satisfied  by 

■n  i  J  i  J 

l 

either  facility  i  or  V  .  In  this  case,  =  |  J\  is  not  guaranteed. 

Remark  2:  Lemma  1  is  only  a  necessary  condition  for  optimality.  For 
sufficiency,  we  also  need  to  ensure,  at  nodes  where  |>Sj  =  |./| ,  that  s  =  0 ,  for  all  j  e  J  , 

i.e.,  there  is  no  unmet  demand.  Consider  the  simple  counter-example  depicted  in  Figure 
1 ,  where  k  =  2.  Suppose  the  incumbent  values  for  bound  coefficients  are  as  shown,  and 
that  c.  =  1  for  all  is  I,  j  e  J  . 

ij  j 
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Facilities 


Customers 


Figure  1.  Bipartite  graph  for  illustrating  Lemma  1. 

This  graph  depicts  an  instance  where  an  optimal  solution  of  x*  =  (1,  1,0)  with  yu  =  20  , 
y2l  =30  and  y22  =  70  does  not  meet  the  invalidity  test  condition  of  |S,|^|,/|,  but  still 
yields  a  suboptimal  solution  to  the  original  MIP. 


An  optimal  solution  to  this  problem  is  x*  =  (1,  1,0),  with  yu  =  20 ,  y21=30  and 
y21  =  70 .  Observe  that  l^l  =  |/| .  However,  this  solution  is  clearly  suboptimal  since  there 
is  unmet  demand  occurring  for  customer  2.  Hence,  a  sufficient  condition  would  also 
demand  that  s  .  =  0 ,  for  all  j  e  .J  ,  whenever  |S|  =  |./| . 

One  approach  to  ensure  optimality  using  only  the  sufficient  condition  is  to  initialize 
bound  coefficients  to  the  average  demand  fulfilled  by  an  opened  facility,  i.e., 

M;.  <—  M\  =  X ,  d  j k  ,  for  all  iel .  In  that  manner,  unmet  demand  will  never  occur. 

We  utilize  Lemma  1  as  another  invalidity  test;  that  is,  at  any  node, 
whenever  |,S'|  ^  |./| ,  we  adjust  the  bound  coefficients  for  all  facilities  that  are  servicing  a 

customer.  The  alternative  invalidity  test  condition  for  Algorithm  3  applied  to  the 
uncapacitated  ^-median  problem  is 

i  e  /;  if  Xy-fy  ^  0  and  |5|  *  M  (3.B.4) 
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We  solve  the  illustrative  example  using  this  alternative  test  condition,  and 
the  following  inputs: 

Input:  Uncapacitated  ^-median  problem  defined  in  previous  section. 

Growth  factor  y  =  (p  =  2 . 

Invalidity  test  condition  |S|  ^  |  j\ . 

M‘  =  YjjdJk  =  540,  Vie/. 

Solving,  Algorithm  3  with  the  alternative  test  condition  yields  the  optimal 
solution  x*  =  (0,  1,  0,  1),  z*=1920. 

C.  COMPUTATIONAL  RESULTS 

We  test  the  algorithms  on  one  randomly  generated  instance  of  the  uncapacitated 
^-median  problem  with  n  -  50  locations,  from  which  k  -  5  must  be  chosen  to  fulfill  the 
demands  of  m  = 1 00  customers.  Data  for  the  model  are  generated  using  the  uniform 
probability  distributions  specified  in  Appendix  A. 

The  test  problem  is  first  solved  using  a  standard  B&B  algorithm,  and  the  optimal 
solution  and  time  taken  to  solve  are  recorded.  We  then  use  the  proposed  algorithms  to 
solve  the  test  problem.  The  measures  of  effectiveness  are: 

1.  The  relative  optimality  gap,  possibly  0%,  between  the  optimal  solution 
derived  from  the  standard  B&B  algorithm  with  baseline  (weak)  coefficients,  and  the 
solution  derived  from  the  proposed  algorithms;  and 

2.  The  time  taken  by  the  proposed  algorithms  to  complete,  as  compared  to 
the  time  taken  by  the  standard  B&B  algorithm  with  baseline  coefficients. 

1.  Solution  by  Standard  Branch  and  Bound  Algorithm 

The  test  problem  is  first  run  using  the  standard  B&B  algorithm,  once  with 
baseline  (weak)  coefficients,  and  once  with  tighter,  empirically  valid,  bound  coefficients. 
Henceforth,  the  test  model  with  weak  bound  coefficients  is  called  the  “base-case.”  Table 
5  shows  the  results. 
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Case 

Bound 

Coefficient,  M j 

Solution 
Time  (s) 

Optimality 
Gap  (%) 

Base-case:  Weak 
Bound  Coefficient 

1,136 

0.0 

Tight  Bound 
Coefficient 

0.3  Yjj  dj 

116 

0.0 

Table  5.  Solution  times  for  standard  B&B  on  the  uncapacitated  ^-median  problem 


As  discussed  in  Chapter  I,  having  weak  bound  coefficients  can  increase  solution 
times.  In  this  instance,  the  time  taken  to  solve  the  problem  using  tight  bound  coefficients 
is  about  10%  of  the  time  taken  to  solve  the  problem  using  baseline  coefficients. 
Naturally,  tight  bound  coefficients  are,  in  general,  unknown. 

2.  Algorithm  1:  Results  and  Analysis 

For  select  values  of  a  and  (3 ,  we  proceed  to  run  Algorithm  1  with  M]  <-  X ,  d. 
Table  6  displays  the  results. 
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Test  Factor 

a 

Decrease 
Factor  (3 

Solution 
Time  (s) 

Number  of  Bound 
Coefficient  Updates 

Optimality 
Gap  (%) 

0.1 

0.1 

2.1 

22,624 

312.2 

2.3 

23,961 

5.1 

1.7 

16,771 

2.5 

6.2 

60,199 

2.7 

0.5 

10.0 

99,719 

2.4 

0.6 

13.3 

97,142 

0.57 

0.7 

36.2 

188,624 

1.2 

0.8 

94.9 

526,975 

0.5 

288.7 

1,071,470 

0.0 

0.2 

0.1 

0.6 

1,551 

315.2 

0.2 

0.3 

714 

57.3 

0.3 

0.3 

1,153 

127.1 

0.4 

1,363 

21.8 

0.4 

1,636 

0.5 

2,642 

31.6 

0.7 

0.9 

6,162 

11.3 

0.8 

3.4 

27,312 

4.7 

0.9 

64.6 

316,758 

1.6 

Base-case  solution  time:  1,136  seconds 

Table  6.  Solution  times  for  Algorithm  1  on  the  uncapacitated  ^-median  problem 


The  algorithm  yields  the  correct  solution  only  when  the  appropriate  parameters 
are  chosen,  in  this  instance,  a-  0.1  and  f3  =  0.9 .  The  solution  time  for  this  setting  of 
parameters  reduces  by  75%  over  the  base-case.  When  a-  0.2,  we  notice  that  the 
optimal  solution  is  never  obtained,  regardless  of  the  decrease  factor.  Since  any 
inefficiency  test  factor  greater  than  0.2  results  in  weaker  testing  criteria,  no  runs  beyond 
a  =  0.2  are  performed. 

For  tests  that  yield  incorrect  solutions,  the  optimality  gap  ranges  from  0.5%  to 
312.2%.  In  certain  cases,  such  as  when  a  =  0.1  and  ft  =  0.3 ,  the  algorithm  takes  1.7 
seconds,  and  yields  a  solution  with  an  optimality  gap  of  2.5%.  That  means  a  reduction  in 
computational  time  of  2.8  orders  of  magnitude  with  respect  to  the  base-case,  with  only 
modest  solution  error. 


27 


We  notice  from  Table  6  that  smaller  values  of  a  and  larger  values  of  (3  tend  to 
yield  correct,  or  close-to-correct,  solutions.  With  a  small  a ,  when  the  inefficiency 
condition  is  met  for  a  given  facility,  the  total  outflow  from  that  facility  is  small  compared 
to  the  value  of  the  controlling  bound  coefficient.  This  indicates  that  the  incumbent  bound 
coefficient  may  be  too  large.  And,  with  a  large  decrease  factor  (3 ,  given  that  potential 
inefficiency  is  detected,  the  bound  coefficient  is  not  reduced  as  much  as  would  be  using  a 
small  (3 .  This  limits  the  chance  of  eliminating  the  optimal  solution  from  the  solution 
space.  Therefore,  using  a  conservative  (small)  test  factor  a  and  a  large  decrease  factor 
/ 3  increases  the  probability  of  Algorithm  1  yielding  the  correct  solution. 

3.  Algorithm  2:  Results  and  Analysis 

Here  we  test  Algorithm  2  for  select  values  of  parameters  8  and  y ,  with  all  M\ 

set  to  the  average  demand  fulfilled  by  an  opened  facility,  i.e.,  Mj  <—  M\  =  X,  d jk  ,  for 
all  i  &  I .  Table  7  presents  results. 


Test 

Factor  8 

Growth 
Factor  y 

Solution 
Time  (s) 

Number  of  Bound 
Coefficient  Updates 

Optimality 
Gap  (%) 

0.03 

1.1 

973.5 

17 

0.0 

1.5 

800.4 

4 

0.0 

0.04 

1.1 

434.1 

4 

0.0 

1.5 

518.8 

17 

0.0 

0.05 

1.1 

291.6 

17 

0.3 

1.5 

434.6 

17 

0.3 

0.10 

1.1 

33.3 

17 

0.3 

1.5 

35.1 

4 

0.3 

0.20 

1.1 

19.1 

11 

2.1 

1.5 

21.1 

4 

1.9 

Base-case  solution  time:  1,136  seconds 


Table  7.  Solution  times  for  Algorithm  2  on  the  uncapacitated  A'-median  problem 
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Algorithm  2  yields  an  incorrect  solution  for  8  >  0.04  .  In  the  instances  when  the 
appropriate  factors  are  chosen  and  the  algorithm  yields  the  correct  solution,  solution 
times  are  reduced  by  up  to  62%. 

For  some  choices  of  parameters  that  do  not  yield  an  optimal  solution,  we  observe 
that  the  solution  times  are  reduced  by  as  much  as  98%,  but  with  an  optimality  gap  of  only 
about  2%. 

One  reason  for  the  considerably  shorter  solution  times  is  that  Algorithm  2  does 
not  restart  the  B&B  search  process  once  the  standard  B&B  algorithm  progresses  beyond 
the  root  node.  Once  Algorithm  2  determines  that  invalidity  test  conditions  are  no  longer 
met  at  the  root  node,  the  bound  coefficients  are  fixed.  However,  it  is  for  precisely  this 
reason  that  Algorithm  2  may  sometimes  yield  suboptimal  solutions.  Notwithstanding  the 
marginally  suboptimal  solutions  that  are  recorded,  we  do  note  the  marked  reduction  in 
solution  time. 

Remark:  The  test  condition  ^  8Mt  implies  that  the  smaller  the  invalidity 

test  factor,  the  easier  it  is  to  meet  the  test  condition  and,  in  turn,  the  more  likely  the 
algorithm  yields  the  optimal  solution.  This  agrees  with  the  test  results,  where  we 
recognize  that  the  algorithm  converges  to  the  optimal  solution  only  for  small  values  of 
invalidity  test  factor  8. 

4.  Algorithm  3:  Results  and  Analysis 

The  computational  run  times  are  recorded  for  Algorithm  3  for  various  values  of 
parameters  8  and  (f),  while  setting  y  =  (p  =  1.1.  Table  8  presents  results. 
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Root  Test 
Factor  8 

Tree  Test 
Factor  (f) 

Solution 
Time  (s) 

Number  of  Bound 
Coefficient  Updates 

Optimality 
Gap  (%) 

0.1 

0.9 

879.5 

43 

0.0 

0.3 

891.7 

52 

0.0 

0.5 

954.8 

49 

0.0 

0.7 

927.5 

49 

0.0 

0.1 

0.7 

710.5 

43 

0.0 

0.3 

834.1 

51 

0.0 

0.5 

856.0 

49 

0.0 

0.7 

923.0 

49 

0.0 

0.1 

0.5 

900.4 

43 

0.0 

0.3 

1,231.2 

50 

0.0 

0.5 

934.7 

49 

0.0 

0.7 

933.1 

49 

0.0 

0.1 

0.3 

892.8 

38 

0.0 

0.3 

1,089.5 

49 

0.0 

0.5 

928.9 

49 

0.0 

0.7 

894.0 

49 

0.0 

Base-case  solution  time:  1,1 

136  seconds 

Table  8.  Solution  times  for  Algorithm  3  on  the  uncapacitated  ^-median  problem 


For  the  various  parameters  used.  Algorithm  3  yields  the  optimal  solution  in  all 
instances.  This  suggests  that  Algorithm  3,  together  with  this  test  condition,  is  insensitive 
to  input  parameters.  Solution  times  are  also  shorter  in  most  cases  tested,  by  up  to  37%. 
This  indicates  that  the  bound  coefficients  have  increased  sufficiently  for  validity,  but  still 
less  than  the  baseline  coefficients. 

As  the  algorithm  tests  for  invalid  bound  coefficients  in  the  B&B  search  tree  and 
restarts  from  the  root  node  whenever  updates  are  made,  we  conjecture  that  Algorithm  3, 
together  with  this  test  condition,  always  yields  the  optimal  solution. 

5.  Algorithm  3:  Results  and  Analysis  with  Alternative  Test  Condition 

We  proceed  to  run  Algorithm  3  using  the  alternative  test  condition  for  selected 
values  for  parameters  y.  Again,  we  set  M;  <—  M\  =  X,  d  Jk  ,  for  all  iel .  Table  9 
displays  the  results. 
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Growth  Factor  y 

Solution 
Time  (s) 

Number  of  Bound 
Coefficient  Updates 

Optimality 
Gap  (%) 

1.1 

5,405.3 

1,446 

0.0 

1.2 

7,208.7 

854 

0.0 

1.3 

5,595.1 

665 

0.0 

1.4 

5,691.2 

665 

0.0 

1.5 

3,879.7 

616 

0.0 

1.6 

6,388.3 

681 

0.0 

1.7 

5,967.8 

665 

0.0 

1.8 

5,877.8 

665 

0.0 

1.9 

3,754.0 

616 

0.0 

2.0 

7,767.9 

692 

0.0 

Base-case  solution  time:  1,136  seconds  ) 

Table  9.  Solution  times  for  Algorithm  3  (alternative  test  condition)  on  the 
uncapacitated  ^-median  problem 


We  notice  that  this  test  condition  yields  the  optimal  solution  for  any  choice  of 
growth  factor  y .  This  is  consistent  with  our  claim  that  bound  coefficients,  at  any  node, 

are  unsuitable  whenever  |*Sj  ^  |  J\ ,  for  instances  where  ctJ  ^  cVj ,  for  any  i*i'  el ,  and  any 

j&J. 

Solution  times  are  excessive  here,  making  this  invalidity  test  condition 
impractical.  Even  with  valid  bound  coefficients,  there  may  be  nodes  yielding  non-integer 
solutions  that  meet  the  invalidity  test  condition  .S’|  ^  |  j\ ,  forcing  the  bound  coefficients  to 
be  increased  and  the  problem  restarted. 
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IV.  TEST  MODEL  II:  SHORTEST-PATH  INTERDICTION 

PROBLEM 


This  chapter  describes  a  SPIP,  states  the  invalidity  test  condition  to  be  used  in 
Algorithm  3  and  solves  a  problem  instance  using  that  algorithm.  As  will  be  explained. 
Algorithms  1  and  2  do  not  apply  here. 

A.  SHORTEST-PATH  INTERDICTION  PROBLEM 

Given  a  directed  graph  G  =  (V,  E )  ,  where  V  defines  a  set  of  nodes  and  E  a  set  of 
directed  edges  connecting  the  nodes,  the  SPIP  involves  a  network  user,  or  defender,  and 
an  attacker.  The  defender  seeks  to  traverse  a  path  of  minimum  length  between  two 
specified  nodes,  s  and  t.  The  attacker  seeks  to  maximize  the  shortest-path  length  of  the 
defender,  and  does  this  by  interdicting  edges  in  the  network  and  increasing  the  effective 
lengths  of  the  interdicted  edges,  or  by  making  those  edges  impassable.  We  assume  that 
interdiction  of  edge  (i,j)  does  not  imply  interdiction  of  edge  (  j,i) ,  if  that  edge  exists. 

Fulkerson  and  Harding  (1977)  study  a  version  of  the  SPIP,  where  they  allow  the 
penalty  on  a  particular  edge  to  increase  as  a  linear  function  of  resources  allocated  to  the 
interdiction  of  that  edge.  Israeli  and  Wood  (2002)  study  the  SPIP  with  binary  interdiction 
effort,  i.e.,  with  fixed  penalties;  they  propose  two  algorithms  to  solve  it.  (See  also  Bayrak 
&  Bailey,  2008.) 

We  solve  the  SPIP  with  binary  interdiction  effort,  formulated  as,  follows: 

Indices 

V  set  of  nodes 

E  set  of  edges 

i  e  V  nodes  in  G,  where  5  is  the  source  node,  t  is  the  sink  node 
k  =  (i,  / )  e  E  edges  in  G,  where  i  is  the  origin  and  j  is  the  destination  node  of  edge  k 
k  g  FS(i )  (k  g  RS(i))  edges  directed  out  of  (into)  node  i 
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Parameters 


ck  cost  to  traverse  edge  k 

Mk  added  delay  if  edge  k  is  interdicted 

Variables 

xk  1  if  edge  k  is  interdicted  by  the  attacker,  0  otherwise 


yk  1  if  edge  k  is  part  of  the  post-interdiction  shortest  path,  0  otherwise 
Formulation 


z*  -  max  min  £  (ck  +  xkMk  )yk 

xgX  y 


s-t-  Z  yk-  S  yk 

keFS(i)  keRS(i) 


1,  i  =  s 

0,  Vi  gV\ {s,t} 
-1,  i  =  t 


y>o 

xg{0,1}|£| 


(4.A.1) 


(4.A.2) 

(4. A. 3) 
(4.A.4) 


Here,  X  defines  the  set  of  resources  available  to  conduct  interdictions. 

The  bound  coefficient  of  interest  here  is  the  penalty  cost  M k .  Setting 

Mk  =  iFlmax  hy!  for  all  k  g  E  ensures  that  interdiction  of  an  edge  makes  that  edge 

k' 

effectively  impassable:  This  is  the  problem  variant  we  wish  to  test. 

The  objective  function  (4.A.1)  seeks  to  maximize  the  minimum  cost  for  the 
defender  to  traverse  the  graph.  Flow-balance  constraints  are  modeled  in  (4.A.2)  to  route 
one  unit  of  flow  from  s  to  t.  While  yk  is  required  to  be  binary,  the  unimodular  property 

of  the  shortest-path  problem  allows  this  integer  requirement  to  be  relaxed  and  replaced 
with  a  non-negativity  requirement. 


34 


One  approach  to  solve  the  SPIP  is  to  take  the  dual  of  the  inner  minimization 
problem  (Fulkerson  &  Harding  1977,  Israeli  &  Wood  2002).  The  following  MIP  results: 


z*  =  max  n.  -  /r  (4. A. 5) 

xeZ.i 

s.t.  TTj  -  ni  <  ck  +  Mkxk,  \/k  =  O',  j)  e  E  (4.A.6) 

ns  ~  0  (4.A.7) 

x  e  {0, 1}|£|  (4. A. 8) 


Here,  nj  denotes  the  unrestricted  dual  variable  associated  with  equation  (4.A.2). 
However,  we  may  assume  ns  -  0  as  the  inner  minimization  problem  has  at  least  one 
redundant  flow-balance  constraint.  And,  once  we  fix  -  0 ,  we  can  also  assume  nj  >  0 
for  all  i  e  V  .  A  standard  B&B  algorithm  will  solve  the  SPIP. 

B.  TEST  CONDITION  FOR  SHORTEST-PATH  INTERDICTION  PROBLEM 

Given  our  interpretation  of  interdiction,  a  value  of  Mk  is  invalid  if  an  edge  is 

interdicted  and,  at  the  same  time,  flow  is  being  sent  on  that  edge.  The  interdiction  plan  is 
obtained  from  the  solution  of  the  MIP  shown  in  (4. A. 5)  -  (4. A. 8)  (i.e.,  the  values  of  the 
variables  xk  ),  and  flows  in  the  network  are  obtained  from  the  dual  variables  of  equation 

(4.A.6)  (which  is  essentially  the  flow  variable  yk  in  the  primal  subproblem).  As  before, 
Vn  is  defined  to  be  the  index  set  of  coefficients  Mi  that  meet  the  invalidity  test  condition 
at  node  n.  For  Algorithm  3,  for  the  SPIP,  the  invalidity  test  condition  is 

k  =  (i,  j )  e  Ivn  if  xk  =  1  and  yk  >  0.  (4.B.  1) 

Increase  factors  y  and  (p  are  defined  so  that  Mk  <—  yMk  for  all  k  g  //’ ,  and  M k  <—  (pMk , 
for  all  k  e  with  n>  1.  Furthermore,  if  either  yMk>\v\max{ck.}  or 

k' 

cpMk  >  \  V\  max{ck,} ,  we  let  Mk  <—  Flmaxfq.J . 
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We  test  the  SPIP  only  with  Algorithm  3.  At  the  time  of  writing,  no  inefficiency 
test  condition  has  been  developed;  hence  Algorithm  1  is  not  implemented.  Also,  it  is 
unlikely  the  root  node  will  yield  integer  interdiction  solutions;  hence  Algorithm  2  is  not 
implemented  as  the  invalidity  test  condition  will  not  be  met  at  the  root  node. 

C.  COMPUTATIONAL  RESULTS 

We  test  Algorithm  3  on  an  instance  of  SPIP  with  |U|  =102  and  |fi|  =  416.  While 

testing  is  limited  to  one  instance  only,  it  should  serve  as  a  proof-of-concept  if  good 
results  are  obtained,  i.e.,  further  investigation  will  be  warranted.  For  the  purpose  of  this 
test  model,  costs  to  traverse  edges  are  randomly  generated  according  to  a  uniform 
distribution  between  1  and  10;  actual  values  used  are  listed  in  Appendix  B. 

1.  Solution  by  Standard  Branch  and  Bound  Algorithm 

The  test  problem  is  first  run  using  the  standard  B&B  algorithm,  once  with 
baseline  (weak)  coefficients,  and  once  with  tighter,  empirically  valid,  bound  coefficients. 
The  results  are  shown  in  Table  12. 


Case 

Bound 

Coefficient,  Mk 

Solution 
Time  (s) 

Optimality 
Gap  (%) 

Base-case:  Weak 
Bound  Coefficient 

f  max{ct,} 

186.1 

0.0 

Tight  Bound 
Coefficient 

0.0l|F|max{cA_,} 

0.4 

0.0 

Table  10.  Solution  times  for  standard  B&B  on  SPIP 
We  again  see  that  having  weak  bound  coefficients  can  be  detrimental  to  solution 

times. 
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2.  Algorithm  3:  Results  and  Analysis 

We  proceed  to  run  Algorithm  3  for  select  values  of  y  and  (p,  where  we  choose 
y  —  (p  as  both  are  increase  factors,  albeit  at  different  stages  of  the  algorithm.  We  also  set 
M[  =  0.001  \V\ max {ck,} ,  for  all  k  eE .  Table  13  presents  results. 


Growth  Factors 

r ,  <p 

Solution 
Time  (s) 

Number  of 
Penalty  Updates 

Optimality 
Gap  (%) 

i.i 

1.1 

38 

0.0 

1.2 

0.6 

20 

3.0 

1.3 

0.4 

14 

3.0 

1.4 

0.5 

9 

3.0 

1.5 

0.3 

9 

3.0 

1.6 

0.5 

13 

3.0 

1.7 

0.3 

11 

3.0 

1.8 

0.4 

11 

0.0 

1.9 

0.5 

10 

3.0 

2.0 

0.4 

8 

3.0 

2.1 

0.5 

9 

3.0 

2.2 

0.4 

8 

3.0 

2.3 

0.5 

9 

0.0 

2.4 

0.4 

8 

3.0 

2.5 

0.5 

10 

0.0 

2.6 

0.5 

8 

3.0 

2.7 

0.4 

8 

0.0 

2.8 

0.6 

12 

3.0 

2.9 

0.6 

11 

0.0 

3.0 

0.3 

6 

3.0 

Base-case  run  time:  186.1  seconds 

Table  1 1 .  Solution  times  for  Algorithm  3  on  SPIP 


Algorithm  3  with  the  described  test  condition  does  not  always  yield  an  optimal 
solution.  However,  for  the  cases  tested,  the  optimality  gap  is  at  most  3%,  and  solution 
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times  are  reduced  by  at  least  99.4%.  While  optimality  is  not  guaranteed  with  this  test 
condition,  the  marked  reduction  in  solution  times  warrants  further  investigation. 

3.  Additional  Remarks 

Xpress-IVE,  with  default  settings  enabled,  solves  the  SPIP  defined  with  baseline 
coefficients  in  1.3  seconds.  This  indicates  that  the  efficiency  provided  by  presolving 
techniques,  cutting-plane  strategies  and  heuristics  cannot  be  ignored.  We  have  not 
carried  out  the  analogous  test  for  our  version  of  the  uncapacitated  k-median  problem; 
however:  some  solvers  actually  transform  the  weak  formulation  into  a  stronger  one 
before  solving,  and  since  we  do  not  know  how  Xpress-IVE  would  handle  our 
formulation,  test  results  would  not  be  meaningful. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  continuous  relaxation  of  a  MIP  containing  standard  big-M  constructs  can  be 
weak  because  large  values  of  M  (“bound  coefficients”)  are  required  to  ensure  the  model’s 
validity.  This,  in  turn,  may  lead  to  slow  convergence  when  solving  the  model  using  a 
standard  B&B  algorithm.  Standard  big-M  constructs  are  used  in  MIPs  such  as  (a)  the 
weak  formulations  of  the  uncapacitated  k-median  problem  and  UFLP,  (b)  the  production¬ 
planning  problem  and  (c)  the  SPIP.  This  thesis  develops  three  algorithms  to  adaptively 
adjust  bound  coefficients  within  the  B&B  tree.  The  goal  is  to  achieve  optimal  solutions 
faster. 

Algorithm  1  starts  with  initially  weak  bound  coefficients  defined  by  standard 
techniques  (baseline  coefficients),  and  reduces  them  within  the  B&B  search  tree. 
Algorithms  2  and  3  start  with  initially  small  and  possibly  invalid  coefficients  and  increase 
them  as  their  respective  B&B  algorithms  proceed.  Algorithm  2  detects  potentially 
invalid  coefficients  at  the  root  node  of  the  B&B  search  tree  only,  while  Algorithm  3 
detects  such  coefficients  throughout  the  B&B  tree.  When  coefficients  are  increased  for 
Algorithms  2  and  3,  the  problem  is  re-solved  from  the  root  node.  Central  to  all  three 
algorithms  is  the  need  for  test  conditions  that  signify  possible  invalidity  or  inefficiency  of 
bound  coefficients.  Algorithms  1,  2  and  3  are  tested  on  an  instance  of  the  uncapacitated 
k-median  problem,  and  Algorithm  3  is  also  tested  on  an  instance  of  the  SPIP. 

The  uncapacitated  k-median  problem  seeks  to  place  production  facilities  so  that 
customer  demands  for  a  single  product  can  be  met  at  minimum  cost,  while  placing  a  limit 
on  the  total  number  of  facilities  to  be  opened.  For  testing  purposes,  we  use  a  weak 
formulation  in  which  each  bound  coefficient  limits  the  total  amount  of  product  that  a 
facility  can  ship  to  customers.  In  the  instances  where  Algorithms  1  and  2  yield  an 
optimal  solution,  we  observe  a  reduction  in  solution  time  of  up  to  75%  and  62%, 
respectively.  Importantly,  when  Algorithms  1  and  2  yield  marginally  suboptimal 
solutions,  substantial  solution-time  improvements  are  recorded.  Algorithm  3  yields  some 
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promising  results:  All  cases  tested  converge  to  an  optimal  solution,  with  up  to  a  37% 
decrease  in  solution  time.  An  alternative  test  condition  for  Algorithm  3  also  yields  an 
optimal  solution  in  all  cases,  but  with  excessive  solution  times. 

The  SPIP  seeks  to  use  limited  resources  to  interdict  (destroy)  arcs  in  a  network  in 
order  to  maximize  the  length  (cost)  of  the  shortest  (least-cost)  path  between  two  nodes. 
In  the  tested  instances,  bound  coefficients  refer  to  the  penalties  required  to  model  edges 
as  effectively  impassable  when  interdicted.  Algorithms  1  and  2  are  not  suited  for  the 
SPIP,  so  only  Algorithm  3  is  tested.  Results  show  an  optimality  gap  of  at  most  3%,  with 
solution  times  reduced  by  at  least  99.4%. 

Notwithstanding  the  promising  results  from  the  test  models,  further  research  must 
be  conducted  before  practical  implementations  can  be  adopted. 

B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

This  research  area  is  open  to  further  developments  of  the  algorithms.  Algorithm  3 
can  be  improved  further  such  that  re-solving  a  problem  is  not  required  every  time  bound 
coefficients  are  updated.  This  will  involve  reexamining  fathomed  B&B  nodes  whenever 
bound  coefficients  are  updated,  in  case  a  node  was  fathomed  incorrectly. 

Another  potential  area  for  future  research  would  combine  the  ideas  behind 
Algorithms  1  and  3  to  create  a  single  algorithm  that  is  able  to  detect  and  update  both 
invalid  and  weak  bound  coefficients.  This  would  allow  the  user  to  initialize  bound 
coefficients  to  nearly  arbitrary  positive  values. 

Other  areas  of  focus  are  fonnulation  of  algorithms  that  are  provably  correct  and 
yield  substantial  reductions  in  run  times.  It  is  possible  this  may  not  include  the  ideas 
used  in  the  adaptive  B&B  algorithms. 

The  algorithms  proposed  in  this  thesis  are  tested  on  instances  of  the  uncapacitated 
/(-median  problem  and  SPIP  with  randomly  generated  data.  Real-world  problems,  such 
as  production  lot-sizing  problems  should  be  identified  and  used  to  test  the  algorithms. 
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APPENDIX  A.  DISTRIBUTION  OF  DATA  FOR  UNCAPACITATED 

^-MEDIAN  TEST  MODEL 

The  uncapacitated  ^-median  problem  tested  comprises  n  =  50  locations,  from 
which  k  =  5  must  be  chosen  to  fulfill  the  demands  of  m  =100  customers.  The  data  used 
are  randomly  generated  according  a  unifonn  distribution  using  random  seed  equal  to  2  in 
Xpress-IVE.  The  details  are  as  follows: 

•  (xFac; ,  yFaci )  denotes  the  (x,y)  coordinates  for  facility  i,  where 

xFaci~  Uniform^ 0,3000)  and  yFaci  -  Uniform^ 0,1000). 

•  (xCustj ,  yCustj )  denotes  the  (x,y)  coordinates  for  customer  j,  where 

xCustj  ~  Unifomi(0, 3000)  and  yCustj  ~  Uniform^ 0,1000) . 

•  ct]  denotes  per  unit  shipping  cost  from  facility  i  to  customer  j,  where 
c(.  j  =  yJ(xFaci  -  xCustj  )  +  (yyFaci  -  yCust  j  j 

•  dj  denotes  units  of  demand  at  customer  j,  where  dj  ~  Uniform^ 100,200) 

•  Pj  denotes  penalty  per  unit  of  unmet  demand  at  customer  j,  where 
P  y  =11  max  {c.} 

FJ 


43 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


44 


APPENDIX  B.  DATA  SET  FOR  SHORTEST-PATH  INTERDICTION 

PROBLEM 


The  SPIP  test  model  uses  an  asymmetric  grid  network  G  =  ( V,E )  with  V  =102 


and  is  =  416 .  Table  14  presents  all  of  the  data  for  the  test  instance,  except  the  facts  that 


s  =  \,t  =  2. 


Arc 

Cost 

Arc 

Cost 

Arc 

Cost 

Arc 

Cost 

(1,3) 

3 

(27,37) 

10 

3 

(72,82) 

1 

(3,13) 

6 

(27,38) 

4 

■EES! 

9 

(72,81) 

3 

(3,14) 

9 

(27,36) 

IBSEI 

3 

(73,74) 

10 

(1,4) 

9 

(28,27) 

10 

ma 

2 

(73,83) 

2 

(4,14) 

10 

(28,29) 

7 

(50,61) 

6 

(73,84) 

1 

(4,15) 

7 

(28,38) 

2 

(50,59) 

6 

(74,73) 

8 

(4,13) 

2 

(28,39) 

4 

(51,50) 

10 

(74,75) 

7 

(1,5) 

6 

(28,37) 

10 

(51,52) 

6 

(74,84) 

4 

(5,15) 

2 

(29,28) 

3 

(51,61) 

9 

(74,85) 

10 

(5,16) 

5 

(29,30) 

4 

(51,62) 

1 

(74,83) 

9 

(5,14) 

9 

(29,39) 

2 

(51,60) 

10 

(75,74) 

6 

(1,6) 

7 

(29,40) 

4 

(52,51) 

5 

(75,76) 

5 

(6,16) 

8 

(29,38) 

1 

(52,62) 

2 

(75,85) 

6 

(6,17) 

9 

(30,29) 

7 

(52,61) 

5 

(75,86) 

7 

(6,15) 

3 

(30,31) 

8 

(53,54) 

8 

(75,84) 

7 

(1,7) 

5 

(30,40) 

(53,63) 

3 

(76,75) 

6 

(7,17) 

3 

(30,41) 

7 

(53,64) 

3 

(76,77) 

9 

(7,18) 

5 

(30,39) 

5 

(54,53) 

6 

(76,86) 

8 

(7,16) 

10 

(31,30) 

2 

(54,55) 

1 

(76,87) 

7 

(1,8) 

10 

(31,32) 

5 

(54,64) 

4 

(76,85) 

10 

(8,18) 

9 

(31,41) 

2 

(54,65) 

7 

(77,76) 

5 

(8,19) 

1 

(31,42) 

6 

(54,63) 

1 

(77,78) 

5 

(8,17) 

8 

(31,40) 

5 

(55,54) 

3 

(77,87) 

9 

(1,9) 

6 

(32,31) 

1 

(55,56) 

3 

(77,88) 

4 

(9,19) 

2 

(32,42) 

4 

(55,65) 

4 

(77,86) 

1 

(9,20) 

6 

(32,41) 

6 

(55,66) 

4 

(78,77) 

1 

(9,18) 

9 

(33,34) 

7 

(55,64) 

6 

(78,79) 

4 

(1,10) 

6 

(33,43) 

9 

(56,55) 

4 

(78,88) 

3 

(10,20) 

9 

(33,44) 

7 

(56,57) 

9 

(78,89) 

8 

45 


Arc 

Cost 

Arc 

Cost 

Arc 

Cost 

Arc 

Cost 

(10,21) 

3 

(34,33) 

10 

(56,66) 

3 

(78,87) 

7 

(10,19) 

8 

(34,35) 

7 

(56,67) 

9 

(79,78) 

4 

(Ul) 

2 

(34,44) 

3 

(56,65) 

7 

(79,80) 

5 

(11,21) 

8 

(34,45) 

9 

(57,56) 

(79,89) 

7 

(11,22) 

6 

(34,43) 

6 

(57,58) 

1 

(79,90) 

1 

(11,20) 

2 

(35,34) 

2 

(57,67) 

6 

(79,88) 

3 

(U2) 

9 

(35,36) 

9 

(57,68) 

7 

EH 

3 

(12,22) 

10 

(35,45) 

8 

(57,66) 

7 

mm 

9 

(12,21) 

1 

(35,46) 

2 

(58,57) 

4 

mm 

6 

(13,14) 

10 

(35,44) 

5 

(58,59) 

3 

(80,91) 

9 

(13,23) 

7 

(36,35) 

9 

(58,68) 

6 

(80,89) 

1 

(13,24) 

1 

(36,37) 

8 

(58,69) 

(81,80) 

2 

(14,13) 

4 

(36,46) 

7 

(58,67) 

8 

(81,82) 

10 

nESBI 

5 

(36,47) 

4 

(59,58) 

1 

(81,91) 

8 

1EEJ 

6 

(36,45) 

1 

(59,60) 

8 

(81,92) 

5 

6 

(37,36) 

1 

(59,69) 

10 

(81,90) 

1 

10 

(37,38) 

3 

(59,70) 

2 

(82,81) 

1 

5 

(37,47) 

10 

(59,68) 

9 

(82,92) 

5 

3 

(37,48) 

9 

10 

(82,91) 

10 

10 

(37,46) 

5 

4 

(83,84) 

7 

6 

(38,37) 

5 

2 

(83,93) 

2 

7 

(38,39) 

8 

1 

(83,94) 

3 

8 

(38,48) 

9 

10 

(84,83) 

10 

4 

(38,49) 

1 

8 

(84,85) 

7 

2 

(38,47) 

7 

3 

(84,94) 

3 

3 

(39,38) 

3 

(61,71) 

3 

(84,95) 

1 

6 

(39,40) 

9 

(61,72) 

3 

(84,93) 

6 

7 

(39,49) 

6 

(61,70) 

(85,84) 

8 

8 

(39,50) 

3 

(62,61) 

5 

(85,86) 

8 

5 

(39,48) 

2 

(62,72) 

6 

(85,95) 

3 

7 

mi 

9 

(62,71) 

8 

(85,96) 

4 

8 

mm 

5 

(63,64) 

9 

(85,94) 

1 

7 

(40,50) 

5 

(63,73) 

5 

(86,85) 

3 

1 

(40,51) 

7 

(63,74) 

1 

(86,87) 

7 

9 

(40,49) 

8 

(64,63) 

5 

(86,96) 

6 

3 

(41,40) 

6 

(64,65) 

(86,97) 

7 

9 

(41,42) 

10 

(64,74) 

3 

(86,95) 

4 

(19,18) 

2 

(41,51) 

2 

(64,75) 

9 

(87,86) 

3 

46 


Arc 

Cost 

(19,20) 

6 

(19,29) 

7 

(19,30) 

3 

(19,28) 

6 

(20,19) 

7 

(20,21) 

10 

(20,30) 

4 

(20,31) 

1 

(20,29) 

5 

(21,20) 

5 

(21,22) 

2 

(21,31) 

9 

(21,32) 

3 

(21,30) 

10 

(22,21) 

1 

(22,32) 

7 

(22,31) 

8 

(23,24) 

4 

(23,33) 

10 

(23,34) 

6 

(24,23) 

5 

(24,25) 

7 

(24,34) 

6 

(24,35) 

2 

(24,33) 

6 

(25,24) 

7 

(25,26) 

10 

(25,35) 

8 

(25,36) 

3 

(25,34) 

2 

(26,25) 

8 

(26,27) 

2 

(26,36) 

3 

(26,37) 

7 

(26,35) 

10 

(27,26) 

4 

(27,28) 

7 

Arc 


1,52 
1,50 
2,41 

.  2,52; 

(42,51) 

(43,44) 


Arc 


64.73 
65,64 
65,66 

65.75 

65.76 

65.74 


66.65 

66.67 

66.76 

66.77 
66,75 

67.66 

67.68 


67.77 

67.78 
67,76 
68,67 
68,69 


68.78 

68.79 
68,77 
69,68 
69,70 

69.79 

69.80 


69,78 

70,69 

70,71 

70.80 

70.81 


70.79 

71.70 
71,72 

71.81 

71.82 

71.80 

72.71 


(90,99) 

(91.90) 
(91,92) 

91.101 

91.102 
91,100 

(92.91) 


Cost 

8 

7 

3 

3 

6 

9 

3 

10 

8 

9 

2 

1 

) 

3 

8 

10 

8 

) 

7 

) 

3 

5 

1 

2 

(92,102) 

2 

(92,101) 

2 

(93,2) 

2 

(94,2) 

3 

(95,2) 

5 

(96,2) 

4 

(97,2) 

10 

(98,2) 

9 

(99,2) 

8 

(100,2) 

10 

(101,2) 

4 

(102,2) 

8 

Table  12.  Data  for  SPIP  test  problem 
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